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Abstract 

Parton Cascade Models (PCM [1-4]), which describe the full time-evolution of a system of 
quarks and gluons using pQCD interactions are ideally suited for the description of jet produc- 
tion, including the emission, evolution and energy-loss of the full parton shower in a hot and 
dense QCD medium. The Landau-Pomeranchuk-Migdal (LPM) effect [5,6], the quantum inter- 
ference of parton wave functions due to repeated scatterings against the background medium, 
is likely the dominant in-medium effect affecting jet suppression. We have implemented a prob- 
abilistic implementation of the LPM effect [7] within the PCM which can be validated against 
previously derived analytical calculations by Baier et al (BDMPS-Z) [8-12]. 
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1. Introduction 

Recent data from the Relativistic Heavy-Ion Collider (RHIC) at Brookhaven have 
provided strong evidence for the existence of a transient Quark-Gluon-Plasma (QGP) 
with the properties of a near ideal fluid and high opacity. The discovery of the QGP has 
been anchored by the measurement of elliptic flow and the suppression of particles with 
high transverse momentum, jet quenching. 

The suppression of high pt hadrons as measured at RHIC [13,14] is a key signature 
of the formation of a QGP and of the modification of jets, collections of collimated 
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high momentum particles, by this matter. The modification of hard probes provides an 
opportunity to learn about the physical properties of the QGP. 

It has been shown that elastic energy loss is insufficient to entirely account for the 
fivefold suppression of high momentum partons as seen in the RHIC data [14]. Medium 
induced radiation must therefore play an important part in this process, furthermore 
at high probe energies quantum effects cannot be neglected [15,16]. The LPM effect is a 
quantum modification to the radiation of highly energetic particles passing through amor- 
phous matter. At high energies forward scattering momentum transfers will be small and 
so the wavefunction of the particle will probe long distances in the material. Interactions 
with the medium can no longer be treated as single scatterings. Instead interactions with 
collections of scattering centers must be treated coherently. This gives rise to interference 
effects which suppress the radiation spectrum. 

In the non-Abelian LPM effect it is the rescattering of bremsstrahlung gluons with 
medium particles which leads to the dominant interference process. During emission 
collinear soft gluons are indistinguishable from their parent partons for some period, 
which can be characterized by a formation time 17 = w/k^_. If this time becomes longer 
than the mean free path, i.e small emission angles and soft radiation, the gluon will 
interact coherently with additional scattering centers. As a result the radiation process 
becomes delocalized, with some characteristic number of scatterers being responsible for 
a single emission. This leads to a suppression of the radiation spectrum and a strong 
length dependence of the energy loss of the emitting parton [8,10]. This effect is believed 
to provide the dominant contribution to the observed suppression of leading hadrons. 

2. The Parton Cascade 

We simulate the evolution of an infinite brick of partonic QGP matter with the 
VNI/BMS parton cascade model [1,2] a fully-relativistic Monte-Carlo transport code 
which numerically solves the Boltzmann equation, 

p^F k (x,p) = J2C i F k (x,p). (1) 

i 

The collision term d includes all possible 2 — » 2 scatterings which may also produce final 
state radiation, 

C i F k (x,p) = ^- j\[dT J \M l \ 2 8 A {P m -P out )D{F k {x 1 p)), (2) 

dTj is the Lorentz invariant phase space for the process j, D is the collision flux factor 
and Si is a process dependent normalization factor. A geometric interpretation of the 
total cross section for each process is used to determine which partons will interact, the 
partons are propagated along classical trajectories. 

The PCM approach allows for the simulation of the evolution of a jet of fast partons 
and the reaction of the medium partons equally. The brick setup allows for a close 
examination of the energy loss processes free from the additional modifications due to 
motion through an expanding and cooling medium. This allows comparison with other 
theoretical models of jet energy loss [19]. 

The medium is initialized with quarks and gluons sampled from thermal distributions 
typically at a temperature of 350MeV, the simulation volume has periodic boundary 
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conditions to simulate an infinite brick of matter. A jet is created by injecting a hard 
parton with energy ~ lOOGeV into the system, this leading parton is then tracked through 
the subsequent scatterings by tagging the hardest outgoing parton from each interaction. 
The strong coupling constant is fixed to a s — 0.3 for scatterings. The final-state radiation 
is produced in an angular ordered parton shower whenever partons emerge from a collision 
with some time-like virtuality [17,18]. 



3. Implementing the LPM effect 



Since the PCM is a Monte-Carlo process it is difficult to naturally include interference 
effects such as the LPM suppression. We implement the LPM effect using a local Monte- 
Carlo routine in the style of Zapp and Wiedemann [7], which reproduces the leading 
BDMPS-Z [8-12] result for light-parton energy loss in a QGP medium AE ~ L 2 . This 
method is particularly appealing since it requires no artificial parametrization of the 
radiative process, it is a purely probabilistic medium induced modification. 

A shower of partons is produced by PYTHIA after an inelastic scattering. The hardest 
radiated gluon is selected to represent the shower as the probe and re-interact with the 
medium. This reflects the dominance of the gluon rescattering in the interference process. 
The formation times 

r°f= E (3) 

branchings 

for each branching during the parton shower leading up to the production of the probe 
gluon are summed. The probe parton is allowed to propagate through the medium and 
rescatter elastically during this time, the remainder of the partons from the radiation 
event propagate spatially but may not interact. Each time the probe gluon rescatters its 
formation time is recalculated as 

T? = ^ (4) 

( k ^ + E"=i cl-m) 

this simulates the emission of the shower from n centers which transfer their momentum 
coherently. After this formation time expires the radiation is considered to have separated 
from the initiating particles and all partons may once again interact and radiate. 



4. Results 



The total energy loss for a hard probe injected into the box is shown as a function of 
time in Fig: 1. This energy loss is only incremented when the probe parton is outside of 
its formation time, when the parton is not in a coherent state with the emitted radiation. 

The quadratic fit in Fig:l reproduces the leading order result from BDMPS [10] 

-AE=^^L 2 , (5) 

8 Ag 

for light partons the time elapsed and distance traveled through the box are intcr- 
changable and we use thermal values for the characteristic momentum transfer fj 2 = 
\g 2 T 2 and the mean free path 1 = 2>TCa<^ s log The elastic energy loss process has 
been found to agree well with theoretical predictions [20]. 
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Fig. 1. The total energy loss for a lOOGeV light quark traversing a QGP medium at T = 0.35GeV, with 
a s = 0.3. The open red circles show the energy loss from radiation with the LPM effect and elastic 
scattering, the closed green circles show the energy loss due to elastic scattering alone. 

We have successfully implemented a local Monte-Carlo version of the LPM effect and 
reproduced to leading order the canonical BDMPS result for the energy loss of the lead 
parton in a jet. The VNI/BMS parton cascade can now be used to examine how energy 
is carried away from the lead parton into the jet and the medium itself. 
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